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PROCEDE DE LOCALISATION D'UN OU DE PLUSIEURS EMETTEURS 

L'invention concerne un precede de localisation d'un ou de plusieurs 6metteurs. 

5 II s'applique notamment a la localisation d'emetteurs au sol a partir 

d'un engin mobile sans avoir de connaissance a priori sur les signaux emis. 

Le domaine technique est notamment celui de la localisation 
passive d'emetteurs. 

10 La figure 1 illustre une localisation aeroportee. L'emetteur est a la 

position (xo,yo,z 0 ). Le porteur a I'instant t k est a la position (x*,y*.zk) et voit 
I'Smetteur sous I'incidence (^/oXo.yo.zo), A(t k ,xo,yo,zo))- Les angles 
^txo,yo,zo) et A(t,x 0 ,yo,z.o) evoluent au cours du temps et dependent de la 
position de l'emetteur ainsi que de la trajectoire du porteur. 

15 Les angles 0(t,xo,yo,zo) et A(t,x 0 ,yo,zo) sont reperes par rapport a 

un reseau de N antennes pouvant etre fixe sous le porteur comme le montre 
la figure 2. 

il existe actuellement plusieurs techniques permettant de 
determiner la position (x m> y m , z m ) d'un emetteur. Ces techniques de 

20 localisation different notamment par les parametres qui sont estimes en 
instantanee au niveau du reseau des capteurs. Elles peuvent etre classees 
de la maniere suivante: 
Utilisation de la goniometrie, 

Ces techniques sont connues et utilisees dans Tart anterieur. 

25 Dans la plupart des cas, elles sont basees sur une goniometrie 1D en 
azimut. Les azimuts 6U=6(f/c,x,„ > y,„ l z m ) associes au m l6me emetteur sont 
mesures pour differents instants f fc En utilisant la position (x k> yk,Zk) du porteur 
a I'instant correspondant k, une position (x mk ,y mk ,z mk ) de l'emetteur m est 
estimee par une intersection sol. La position (x k ,yi<,Zk) du porteur est donnee 

30 par un GPS, son orientation est obtenue par un compas dans le cas d'un 
porteur terrestre et par une centrale de navigation dans le cas d'un aeronef. 
A partir de toutes les positions (Xm^y^Zm^ la methode effectue une 
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extraction de donnees permettant de determiner les M positions dominantes 
(Xm,ym,Z/n) des emetteurs incidents. La localisation est obtenue par 
triangulation ou par intersection sol (goniometrie 2D). L'inconvenient des 
techniques de triangulation est qu'elles necessitent un defilement important. 
5 D'autre part, les techniques de goniometrie doivent utiliser un reseau de 
capteurs non ambigue pour fournir les incidences. Ceci a pour inconvenient 
de necessiter une table de calibration et de limiter la taille du reseau de 
capteurs et par consequent de fournir des incidences limitees en precision. 
Utilisation de la difference de phase entre 2 capteurs eloignes, 
10 La difference de phase Aq^t^xo^o^o) entre capteurs depend de la 

positions des 2 capteurs ainsi que de Pincidence (6(^Xo,yo,Zo), Ait^Xo^o^o)) 
de Pemetteur. Cette phase qui depend du temps est directement liee a la 
position (xo,yo,zo) de Pemetteur. En consequence, en etudiant la fonction du 
temps A<p(t,x 0 ,yo,z.o) il est possible d'en deduire la position (xo,y 0 ,zo) de 
15 Pemetteur. Dans cette famille d'application les 2 capteurs sont eloignes pour 
augmenter la precision de la mesure de la phase. Ceci a pour inconvenient 
de faire varier la difference de phase A<p(t,Xo,yo,zo) en fonction du temps sur 
plus de 2n et la technique necessite alors une etape permettant de derouler 
la phase sur plus de 2n. D'autre part dans cette technique la phase est 
20 mesuree en effectuant directement une intercorrelation entre 2 capteurs ce 
qui ne permet pas de traiter le cas multi-emetteurs. 
Utilisation de la mesure de la frequence porteuse de I'emetteur, 

Ces techniques exploitent le fait que la frequence porteuse 
estimee est la somme de la frequence porteuse de I'emetteur et du decalage 
25 doppler du a la Vitesse de deplacement du porteur. Le d6calage doppler a 
Tavantage de dependre de la position (x 0 ,yo,zo) de Pemetteur et d'etre aussi 
une fonction du temps Af(f,x 0 ,yo,Zo). En consequence en etudiant la fonction 
du temps Af(txo,yo,z 0 ) il est possible d'en deduire la position (xo,yo,z 0 ) de 
Pemetteur. La mesure de ce decalage doppler presente toutefois comme 
inconvenient de necessiter des emetteurs ayant des formes d'ondes 
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particulieres. Cette mesure de frequence peut se faire par des techniques 
cycliques supposant que le signal emis est non circulaire. 
Utilisation des temps de propagation, 

Ces techniques exploitent les differences de temps de propagation 
5 entre aeriens (TDOA ou Time difference of arrival) qui sont directement liees 
aux distances respectives de I'emetteur aux differents aeriens et done a la 
position (xo,yo,2o) de I'emetteur. En utilisant au moins trois aeriens 
suffisamment espaces, il est possible de deduire ia position (xo,yo,zo) de 
Temetteur par localisation hyperbolique. L'inconvenient de ces techniques est 
10 qu'elles ne peuvent etre mises en oeuvre en contexte mono porteur en raison 
des espacements considerables requis entre aeriens. D'autre part dans ces 
techniques, la difference de temps est mesuree en effectuant directement 
une intercorrelation entre 2 capteurs, ce qui ne permet pas de traiter le cas 
multi-emetteurs. 

15 

Le procede selon I'invention repose notamment sur une 
nouvelle approche d'estimation directe des positions (x m , y m , z m ) de chacun 
des emetteurs a partir d'une analyse parametrique du signal multi-voies a 
divers instants tk sur une duree Dt. L'analyse parametrique a notamment 
20 pour fonction supplemental de separer les differents emetteurs a chaque 
instant tk. On associe ensuite les parametres d'un meme emetteur issus des 
differents instants tk pour finalement localiser chacun des emetteurs. 

Uinvention conceme un procede de localisation d'une ou de 
25 plusieurs sources, la ou lesdites sources etant en mouvement par rapport a 
un reseau de capteurs, le procede comportant une etape de separation des 
sources afin d'identifier les vecteurs directeurs associes a la reponse des 
capteurs a une source d'incidence donnee. II est caracterise en ce qu'il 
comporte au moins les etapes suivantes : 

• associer les vecteurs directeurs a^.-.a^ du m i6me emetteur obtenus 
respectivement aux instants t r ..t K , 
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• localiser le m i6me Smetteur a partir des vecteurs a 1m ...a Km associes. 

Le procede selon Pinvention presente notamment les avantages 

suivants : 

5 • il permet de localiser en plus de la position en (x, y, z) d'un emetteur son 
vecteur vitesse, 

• il s'applique lorsque Ton est en presence de un ou plusieurs 6rnetteurs 
incidents, 

• sa mise en oeuvre ne necessite pas de connaissances particulieres sur le 
10 signal emis, 

• il permet d'utiliser un r§seau de capteurs ambigue c'est-a-dire plusieurs 
incidences sont associees a la mdme r^ponse du reseau qui ont 
Pavantage d'etre grand et ainsi d'etre plus robuste aux phenomenes de 
couplage entre aeriens ou plus generalement aux erreurs de moderation 

15 du reseau d'aeriens, 

• il peut etre mis en oeuvre sur des reseaux calibres en (9,A). 

• il peut etre mis en oeuvre sur des reseaux a antennes a diversites 
d'amplitude comme les antennes colocalisees : reseau avec des dipoles 
de meme centre de phase et ayant des orientations differentes. 

20 

D'autres caracteristiques et avantages de Tobjet de la presente 
invention apparaTtront mieux a la lecture de la description qui suit donnee a 
titre illustratif et nullement limitatif a la lecture des figures annexees qui 
representent : 

25 • la figure 1 le schema de principe de la localisation d'un emetteur de 
position au sol au moyen d'un aeronef, 

• la figure 2 la relation entre un reseau d'antennes et Pincidence d'un 
emetteur, 

• la figure 3 un schema general expliquant le fonctionnement du procede 
30 selon Pinvention, 
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• les figures 4, 5 et 6 des exemples de mise en oeuvre du precede selon 
invention. 

Afin de mieux faire comprendre I'objet de la presente invention, la 
5 description qui suit est donnee £ titre illustratif et nullement limitatif pour 
localiser plusieurs §metteurs disposes au sol au moyen d'un reseau de 
capteurs equipant un aeronef en mouvement. Un tel systeme est par 
exemple decrit a la figure 1. L'aeronef est equipe d'un processeur adapts a 
mettre en oeuvre les etapes du procede selon Pinvention. 
10 Le procede peut aussi etre mis en oeuvre dans le cadre de 

vehicule en mouvement au sol. 

La figure 3 repr6sente, dans un diagramme temps-amplitude du 
signal, le signal x(t) compose d'une combinaison des signaux des emetteurs 
a differents instants t-i, t2, ...ti<. Sur cette figure les differentes etapes mises 
15 en oeuvre 3 savoir la separation des emetteurs SE et I'estimation 
parametrique EP, I'association des parametres de chaque §metteur, la 
localisation d'un emetteur sont resumees. 

En presence de M emetteurs, le procede dispose, & Tinstant t en 
sortie des N capteurs du reseau, du vecteur x(t) representatif du melange 
20 des signaux des M emetteurs. Autour de I'instant t k , le vecteur x(f+fc) de 
dimension Mx1, representant le melange des signaux des M emetteurs, 
s'exprime de la maniere suivante : 

M ( 

x(f+fc) = £ aCOfon, Akm) s m (t+t k )+ b(Ht k )=A k s(t+t k ) + b(t+t k ) pour \t\<At/2 v 

ou b(f) est le vecteur bruit suppose gaussien, a(9,A) est la reponse du reseau 
de capteurs a une source d'incidence (0,A), A*=[a(0/ci, A/ri)..- a(BkM, A^)], 
25 s(f)=[si(0...SA^0] T i Qkm=0itfoX mi y mi z m ) et A km =A(tfoX m9 y mi z m ). Dans ce modele 
la matrice A* de melange depend de I'instant fod'observation. 
Le vecteur directeur de Tincidence correspondant au nri* me emetteur a 
I'instant t k 
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aw* a(Q km , Akm)=^(t h x m ,y m3 z m ) du m ,eme emetteur 



(2) 



est une fonction connue de t k et de la position de I'emetteur (x m ,y m ,z m ). 



Le procede selon 1'invention comporte au moins les etapes 



suivantes : 

5 1. estimer un ou plusieurs parametres associes a la position de la source, 
par exemple les vecteurs directeurs, les incidences, la position, etc. et 
separer les M emetteurs pour les differents instants t^, ce qui consiste a 
identifier les vecteurs directeurs d'incidence a km pour (1<m<M). Cette 
premiere etape est par exemple effectuee par des techniques de 
10 separation de sources connues de I'Homme du metier, 

2. associer les parametres estim6s pour le m^ me emetteur, par exemple en 

associant les differents vecteurs directeurs d'incidences, a-| m a|< m 

obtenus respectivement aux instants t^ ,...tk, 

3. localiser le m idme emetteur & partir des vecteurs associes. 
15 Etape dissociation 



procede possede a Tinstant t k les M signatures a/cm pour (1 <m<M). A I'instant 
la separation de source donne les M vecteurs b,- pour (1 <i<M). L'objectif 
de ce suivi est de determiner pour le m^ me emetteur, I'indice i(m) qui 
20 minimise Tecart entre a*™ et bi( m ). Dans ce cas on en deduira que a k+ i, m = 
b/Ym;. Pour effectuer cette association on definit par exemple la distance entre 
deux vecteurs u et v par : 



En presence de M emetteurs et apres separation de sources, le 




(3) 



Ou u H est le transpose conjugue du vecteur u. 
Dans ces conditions Tindice i(m) verifie : 

dfa/cm, bj (m) ) = min [d(a^ &/)] 



(4) 
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done 



d(a km , b i(m )) = min[1- 




Dans cette association on considere une fonction p m associee au m l{ 



5 Au fil de Tassociation on obtient pour chaque emetteur m et pour 1<m<M, la 
fonction p m (t). Cette fonction a notamment pour objectif d'eliminer les instants 
t k dont la valeur p m (t k ) paraTt trap eloignSe d'une interpolation de la fonction 
p m (f), e'est-a-dire que Ton elimine les instants aberrants qui peuvent etre 
associes a d'autres emetteurs. On definit une zone de tolerance +/-A autour 

10 de la courbe definie par la fonction p m (f/c). Cette zone de tolerance dependra 
de la precision d'estimation des vecteurs directeurs a*™. En particulier en 
presence de M=1 source la zone sera de I'ordre de A=3/V^At (ou At est le 
temps elementaire d'estimation parametrique illustre fig.3 et B est la bande 
instantanee du signal x(t) ). 

15 Les etapes de cette association pour K instants t k sont par 

exemple les suivantes : 

Etape ASE - 1 : Initialisation du processus a k=2. Le nombre M d'emetteurs 
initial est par exemple determine par un test de detection du nombre de 
sources a Tinstant to connu de THomme du metier, 
20 Etape ASE - 2 : Pour 1 <m<M determination des indices i(m) en appliquant 
Tequation (4) et en utilisant le vecteur a ktm avec 1<A77<M et les vecteurs b/ 
identifies a Tinstant f/c+i pour (1</<A^), 

Etape ASE - 3 : Pour 1<m<M effectuer Toperation ak+im= b/^; f 
Etape ASE - 4 : Incrementation /c<-/c+1 et si k<K retour a Tetape ASE-1 , 
25 Etape ASE - 5 : A partir de la famille d'instants 0>={ U<...< t K }, eliminer les / 
instants f/eO tel que les coefficients p m (f/) n'appartiennent pas a la zone 
delimitee par la courbe d'interpolation des p m (tk) et la zone de tolerance A. 



emetteur: 



Pm(fo) - d(a/OT7, Bom) 



(5) 
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On eliminera aussi les instants t k ou |p m (t k )- p m (t k -i)|<A. Apr6s ce tri la 
nouvelle famille d'instants est <!>={ U<...< t t } et on pose K=l. 

A la fin de ces etapes, le proc6de a determine les vecteurs ai m 

a/<h,associ§s au ni §me emetteur. 

5 Localisation d'un emetteur 

Le proced§ determine la position du rrt 6me emetteur a partir des 
composantes des vecteurs ai m jusqu'a a*™. Ces vecteurs a*™ ont la 
particularity de dependre de I'instant t k et surtout de la position (x m ,y m ,z m ) de 
Pemetteur. En particulier pour un reseau compost de N=2 capteurs espaces 
10 d'une distance de d dans I'axe du porteur le vecteur verifie akm: 

1 

= ^p^|cos(^ t5 x OT ,7 n>9 zJ)cos(A(/ i ,x m , > ; m ,z OT )) 

La valeur 1 de la premiere composante correspond au capteur de 
reference. D'apres la figure 1, I'incidence (^.x^.y^Zm), A^Xm^Zn,) ) peut 
etre directement calculee a partir de la position (Xk,y k ,z k ) du porteur a I'instant 
t k et la position (Xm,y m ,z*n) de I'emetteur. 

15 Etape de transformation du vecteur 

Selon une premiere variante de realisation, le procede comporte 
une 6tape de correction des a*™, la mesure des vecteurs directeurs a^n est 
gen6ralement obtenue a un facteur complexe pres indetermine. Selon cette 
premiere variante, le procede comporte une etape qui consiste a changer la 

20 reference de phase du vecteur directeur mesure en se ramenant au 
barycentre de phase (defini a un coefficient scalaire constant pres que Ton 
peut fixer a 1). Cette operation est realisee, par exemple, en estimant le 
coefficient de correction determine par la transformation suivante des a/on en 
aW 



= a(f fo x m ,y /n ,z m ) 



(6) 
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Le coefficient de correction n'est pas totalement determine par 
cette expression compte tenu de redetermination d'ordre N de la racine 
complexe. Un suivi de revolution de phase pendant la periode d'observation 
5 est done effectu6. 

Le coefficient complexe etant defini a un facteur pres parmi les N 
racines N l6me de I'unite, le suivi de phase consiste a fixer arbitrairement le 
premier coefficient de correction (en prenant la racine 1 par exemple), puis a 
determiner a chaque nouvelle iteration k+1, le coefficient qui minimise les 
10 ecarts de phases moyens entre le vecteur directeur recentre a k+1 et le 
vecteur recentre & I'instant k. 

Le critere de minimisation, pour des mesures a la m§me 
frequence, peut etre egal a : 



15 

oti les ak+i sont les vecteurs directeur recentres avec le coefficient de 
correction determine arbitrairement par Tune quelconque des racines Nemes 
de I'expression. Pour des mesure a des frequences differences, il est 
possible de comparer les phases des composantes des deux vecteurs 

20 directeurs en les corrigeant d'une puissance donnee par le rapport de ces 
deux frequences. 

Si Ton considere les vecteurs b^m = a'km , il est alors possible de 
comparer cette mesure a la valeur theorique b(tk i Xm,ym 1 Zm) pour laquelle le 
vecteur directeur theorique a^Xm.ym.Zm) est calcule pour une origine 

25 consideree au barycentre (geometrique) theorique de phase (lieu 
geometrique pour lequel la somme theorique des differences de phase 



£x\2& - mod( 



arg( 



(8> 
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s'annule). Ce lieu ne coincide pas (en general) avec le centre de phase du 
reseau. 

Selon une autre variante de realisation, le procede comporte une 
etape de transformation du vecteur en un vecteur b km dont les 
5 composantes sont formees a partir des composantes du vecteur a/cm ■ En 
particulier, le procede constant par exemple le vecteur de dimension (N- 
1)x1 en choisissant un capteur de reference en n=i : 

(9) 

= b(t k ,x m ,y m ,z m ) 
ou a km (i) est la / 6me composante de 

Les composantes de b km correspondent dans ce cas aux rapports des 
composantes du vecteur et du vecteur a^i). 
10 Ainsi dans I'exemple de liquation (6) en fixant /=1 on obtient : 

b km = |^exp^^ (10) 
= a/cm(2)/ a/cm(1) 

Sachant que les vecteurs directeurs a km sont estimes avec une 
certaine erreur e/™ tel que a(f/ 6 x m ,y m> 2 m )+ e^, on peut en deduire qu'il 
en est de meme pour le vecteur transform^ b^ de (9) . 
Etape de maximisation d'un critere de correlation 
15 Sachant que le vecteur a*™ est une fonction de la position 

(x m ,ym,Zm) de I'emetteur il en est de meme pour le vecteur b^. Le procede 
comporte une etape de maximisation d'un critere de correlation vectorielle 
normalise L^x.y.z) dans I'espace (x,y,z) de position d'un emetteur ou 
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U(x,y,z)= 7 — 5 — K, — - — 7 



(11) 



Avec 



'Im 



_ Km _ 



= v K (Xm,y m ,2m) + w K , v K (x,y,z) = 



b(/ 13 x,y,z) 
b(t K9 x 9 y 9 z) 



et wk = 



w 



lm 



w 



Km 



Le vecteur bruit w* a pour matrice de covariance R=E[wk wk h ]. En 
faisant I'hypoth^se que la matrice R est connue, le critere peut etre envisage 
avec une technique de blanchiment. Dans ces conditions on obtient le critere 
L.K(x,y,z) suivant : 

(12) 



L K '(x,y,z)= 



\b K H R 1 y K (x>y>z)\ 



Avec R=E[wk Wk h ] 

5 • II faut remarquer que Ies criteres des equations (11) et (12) sont egaux 
lorsque R=a 2 l, c'est a dire lorsque Ies erreurs sont considerees de niveau 
egal sur tous Ies capteurs et independantes entre capteurs. Le critere 
LK(x,y,z) de I'equation (12) est done valable pour un bruit de modele w K de 
statistiques blanches. 
10 • Les criteres LK(x,y,z) et l_K(x,y,z) sont compris entre 0 et 1 et verifient 
LK(x f y f z)= L/c'(x,y,z)=1 pour la position (Xm,y m ,z m ) du n? l6me emetteur. Cette 
normalisation permet de fixer un seuil de bonne localisation r\. Ainsi tous les 
maximums (x m ,y m ,z m ) de Ufoy.z) qui verifient UfXm.ym.Zm^ *1 sont 
consid§res comme des bonnes localisations. Le seuil peut etre fixe en 
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fonction d'une connaissance approchee des statistiques de w/c. 

• Les criteres L.K(x,y,z) et LK(x,y,z) ont I'avantage de pouvoir mettre en 
oeuvre une technique de localisation en presence d'un reseau de capteurs 
calibres dans I'espace (8,A). Sachant qu'a I'instant t k on connait la relation 

5 analytique liant I'incidence (e(fr,x,y,z), A(f/c,x,y,z)) de I'emetteur a sa position 
(x,y,z), on peut alors deduire a partir de Tincidence (8(f/c,x,y,z), A(fo,x,y,z)) le 
vecteur a(f/oX mi y m ,z m )= a(e(fc,x,y,z), A(f/f,x J y,z)) en realisant une interpolation 
de la table de calibration (relative aux antennes calibrees). Remarquons 
cependant que ce precede est insensible a un biais en phase (en raison du 
10 critere de correlation vectorielle). 

• Ces criteres permettent aussi de tenir compte de la phase et de I'amplitude 
des composantes de a(9,A). La methode peut done etre envisagee avec des 
reseaux a antennes colocalisees a diversite de diagramme. 

II faut remarquer que dans un contexte aeroporte la connaissance de 
15 I'altitude h de Tavion permet de reduire le calcul du critere dans I'espace de 
recherche (x,y) en posant z=/?. Dans I'exemple des equations (6) et (10) le 
vecteur v/c(x,y,z) s'ecrit de la maniere suivante : 

(13) 



v K (x,y,z) = 



exp^y 2xjcos(0(t x ,x,y, z)) cos(Afo ,x, y, z))j 
exp|y 2tt ^ cos(0(t K , x, y, z)) cos(A(^ , x, y, z))j 



Dans ce procede il est possible d'envisager d'initialiser 
ralgorithme a K=Kq puis ensuite de calculer de fagon recursive le critere 
20 Ufoy.z). Dans ces conditions L.K(x,y,z) se calcule recursivement de la fagon 
suivante : 



U+i(x l y,z)= 



(x 9 y 9 zf 



(14) 



ou 
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a K+ i(x,y,z)= a/<(x,y,z) + b K +i m H b(f K+ i,x,y,z) 
y K+ i(x,y > z)=Y K (x J y > z) + b(f K+ i,x,y,z) H b(f K+ i,x,y,z) 

Les coefficients a K+ i(x,y,z)= a^x.y.z) yK + i(x,y,z)= y^x.y.z), Pk+i=Pk 
sont des spectres intermediates permettant de calculer L-K+i(x,y,z). 

Lorsque les vecteurs b(fK+i,x,y,z) et bi< m sont de normes 
constantes egale a p la relation de recurrence de Tequation (14) devient : 

i / x, kW*»y»*)T (15) 

u+l(x>y ' z) - PiK+lf 

ou 

a K+ i(x,y,z)= aK(x,y,z) + b K +i m H b(fo + i,x,y,z) 

5 Le proced6 est decrit jusqu'ici en supposant que les emetteurs ont 

des positions fixes. II peut facilement s'etendre au cas de cibles mobiles de 
vecteur vitesse (Vxm.Vym.Vzm) pour lesquelles on dispose d'un modele 
devolution. Dans ces conditions I'incidence du m® me emetteur se 
parametrise de la maniere suivante : 

fc) (16) 

et Aiari=A( tk, X m -V xm tk ,y/n-Vym t k , Znr^zm tk) 

10 ou (x m ,y m ,z m ) est la position de I'emetteur a I'instant to et (Vxm.Vym.Vzm) les 
composantes de la vitesse de I'emetteur a I'instant to. Dans ces conditions le 
vecteur b km de Tequation (9) est parametre par (x^.y^.z^) et (Vxm.Vym.Vzm) de 
la maniere suivante : 

bfcm = b((foX/n,ynj|Z m »VxmiVyniiVznT» )+WAm (17) 

De maniere naturelle les criteres de localisation Lk et Lk des 
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equations (11) et (12) ne sont plus parametres seulement par (x,y,z) mais 
aussi par (v x ,v y ,v 2 ). Le procede consiste done a maximiser le critere L/^x.y.z, 
v x ,v y ,v 2 ) en fonctions des 6 parametres (x,y,z, v Xj v y ,v z ). 

Le procede peut s'appliquer a un tres grand nombre de mesures. 
5 Dans ce cas, le proced§ comporte une etape de reduction de la complexite 
numerique de calcul (qui est fonction du nombre de mesures) en diminuant 
K. Le procede prevoit d'effectuer sur les mesures elementaires les 
traitements suivants : 

- decimation des instants fo, en eliminant les instants voisins pour lesquels 
10 revolution de la courbe p m (fc) n'est pas significative, 

filtrage (lissage des mesures qui sont les vecteurs directeurs) et sous 
echantillonnage, 

- les mesures sont ensuite fusionnees sur une duree definie (extraction par 
association de vecteur directeur pour produire une mesure de synthese). 

15 Recapitulation des etapes du procede 

Le proc6de de localisation de plusieurs emetteurs utilisant K instants k 
peut se resumer par les etapes suivantes : 

Etape n°1 : Identification des vecteurs a^n pour (1 <m<M) aux K instants fk en 
appliquant par exemple une technique de separation de sources et 
d'identification de sources comme decrit dans les references [2] [3]. 
Etape n°2 : Association des vecteurs a^ m jusqu'a BKm obtenus aux instants 
respectifs fi... f K associe au m^ me emetteur pour ]<m<M en appliquant les 
etapes ASE-1 jusqu'a ASE-5 decrit ci-dessus. 
Etape n°3 : Initialisation du processus a /77=1 

Etape n°4 : Transformation des K vecteurs akm en des vecteurs b/<m comme 
le suggere Tequation (9). 

Etape n°5 : Calcul et maximisation du critere LK(x,y,z) de Pequation (11) pour 
obtenir la position (x m ,y m ,Zm) du rri eme emetteur. 
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Etape n°6 : Incrementation m<— et si m<M retour a I'etape n°3 
Afin d'affiner I'estimation de la position (x m ,y m ,y m ) des emetteurs les etapes 
du procede peuvent etre realisees de fagon iterative de la maniere suivante : 
Etape n°7 : Identification des vecteurs b, pour (1</<M) a ('instants f K +i en 
5 appliquant par exemple une technique de separation et d'identification de 
sources comme decrit dans les references [2] [3]. 

Etape n°8 : Pour *\<m<M d6termination des indices i(m) en appliquant 
I'equation (4) et en utilisant le vecteur a^n et les vecteurs b/ pour (1<i<M). 
Etape n°9 : Pour 1 <m<M est effectue I'operation a^+i m = b/^; 
10 Etape n°10 : Pour 1<m<M calcul du critere LK+i(x,y,z) de fa?on iterative en 
utilisant les equations (14) et (15) et minimisation de LK+i(x,y,z) pour obtenir 
la position (Xm.ym.Zm) du m l * me 6metteur. 

Etape n°11 : Si Ton decide de continuer pour etre plus precis et moins 
ambigue le procede retourne a T6tape n°7. 

15 

Exemple de mise en ceuvre du procede 

Les simulations ont ete r§alisees avec un reseau de A/=2 capteurs 
alignes dans Taxe du porteur avec d/X=3. Comme d/A.=3 une methode 
effectuant une goniometrie aux instants serait completement ambigue et ne 

20 permettrait pas de faire par la suite des triangulations pour effectuer la 
localisation de Temetteur. Sur les figures 5, 6 et 7 correspondant aux criteres 
de localisation pour K=3, 7 et 16 est trace le pseudo-spectre Z-K(x,y) a 
maximiser permettant de determiner la position de Temetteur dans Tespace 
(x,y). Sachant que si Temetteur se situe en (xo.yo) alors L| C (xo J yo)=1 I on en 

25 deduit que les courbes iso-niveaux /.K(x,y)=0.99 caracterisent la largeur du 
lobe principale. Remarquant que la precision de localisation depend de la 
largeur de ce lobe, on en deduit d'apres les figures 5, 6 et 7 que plus K est 
important et meilleur sera la precision de localisation. 
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REVENDICATIONS 

1 - Procede de localisation d'un ou de plusieurs sources, ia ou lesdites 
5 sources etant en mouvement par rapport a un reseau de capteurs, le 

proc6de comportant une etape de separation des sources afin d'identifier les 
vecteurs directeurs associes a la reponse des capteurs a une source 
d'incidence donnee, caracterise en ce qu'il comporte au moins les etapes 
suivantes : 

10 • associer les vecteurs directeurs a^-a^ obtenus pour le rr# me 
emetteur et respectivement pour les instants t r ..t K , 
• localiser le m i6me emetteur a partir des vecteurs a 1m ...a Km associes. 

2 - Procede selon la revendication 1 caracterise en ce que I'etape 
15 ^association comporte au moins les etapes suivantes : 

Etape ASE - 1 : initialiser le processus & k=2, 

Etape ASE - 2 : pour 1<m<M determiner les indices i(m) en utilisant la 
relation d(a/ OT7 , bi (m) ) = mm^a*™, to,)], le vecteur a k , m et les vecteurs b/ 

identifies a Tinstant t k ^ pour (1 <i<M), etablir une fonction p m (tk) = d(a km , a om ) 
20 Etape ASE - 3 : pour 1 <m<M effectuer Toperation a k +i m = bi (m) , 

Etape ASE - 4 : incrementer /c<-/c+1 et si k<K retourner a l'6tape ASE-1 , 
Etape ASE - 5 : a partir de la famille d'instants 4>={ ^<...< t K }, ainsi 
obtenue, extraire les instants f/ qui n'appartiennent pas a une zone definie 
par la courbe p m (tk) et une zone de tolerance. 

25 

3 - Procede selon la revendication 1 caracterise en ce que Tetape de 
localisation comporte au moins les etapes suivantes : 

maximiser un critere de correlation vectorielle normalise L k (x, y, z) dans 
Pespace (x, y, z) de position d'un emetteur avec 
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LK(x,y,z)= 



Pjc H v K (x,y,z)\ 



Avec 



b 



Im 



h(t x ,x,y,z) 



= v K (x m ,y m ,z m ) + w K , v K (x,y,z) = 



b 



w 



Im 



et wk = 



w 



5 0C1 w/f est le vecteur bruit pour toutes les positions (x, y, z) d'un emetteur. 

4 - Procede selon la revendication 3 caracterise en ce que le vecteur b K 
comporte un vecteur representatif du bruit dont les composantes sont 
fonctions des composantes des vecteurs ai m ... aKm- 



5 - Procede selon la revendication 3 caracterise en ce qu'il comporte une 
etape ou Ton determine la matrice de covariance R=E[w K wk h ] du vecteur 
bruit et en ce que Ton maximise le critere 



15 6 - Procede selon la revendication 5 caracterise en ce que revaluation du 
critere L^x.y.z) et/ou du critere L K ' (x.y.z) est recursive. 

7 - Procede selon Tune des revendications 1 a 6 caracterise en ce qu'il 
comporte une etape de comparaison des maximums avec une valeur seuil. 

20 



10 




Ib/R- 1 v K (*,y,4 2 
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8 - Procede selon Tune des revendications 1 a 7 caracterise en ce que la 
valeur de Kest initialement fixee a Kq. 

9 - Procede selon Tune des revendications 1 a 8 caracterise en ce que les 
5 emetteurs a localiser sont mobiles et en ce que le vecteur considere est 

parametre par la position de I'emetteur a localiser et le vecteur vitesse. 



10 
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